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ABSTRACT 

Modem (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high- 
redshift galaxies (z > 5). However, at these redshifts the cosmic microwave background (CMB) temperature 
is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the 
local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of 
high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission 
(e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; 
and (ii) it is a non-negligible background against which the line and continuum emission are measured. We 
show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift 
galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to 
compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function 
of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift 
galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely 
underestimated for cold systems if the impact of the CMB is not properly taken into account. 
Subject headings: galaxies; ISM - galaxies: evolution - sub-millimeter: galaxies, ISM. 



1. INTRODUCTION 

Modern (sub-)millimeter interferometers allow us to rou- 
tinely measure the gas and dust content of very high redshift 
galaxies (z > 5), giving precious insight into the star forma- 
tion properties and physical state of the interstellar medium 
(ISM) in early galaxies. With ALMA, for example, thanks 
to the combination of increased sensitivities and the nega- 
tive ^-correction in the (sub-)millimeter, it will be possible, 
for the first time, to detect the dust continuum and CO lines 
from galaxies with total luminosities that are close to that of 
the Milky Way. Instead of detecting only high-luminosity 
starbu rsts at high red shifts (i.e. sub-millimeter galaxies; 
e.g. Blain et alj 12002 '). we will gain access to the larger 
population of low-luminosity (with infrared luminosities < 
10'^ Lq) galaxies. In the local Universe, the bulk of dust 
in normal star-forming galaxies of moderate in frared lumi- 
nosity typically has temperatures of ~ 20 K (e.g. lSmith et alJ 
120121) . However, at higher redshifts, the cosmic microwave 
background (CMB) temperature approaches and can even sur- 
pass this temperature. This means that we must consider the 
effects of the CMB on our observations. 

The temperature of the CMB at any redshift z is given by: 

TcMBiz) = T^MB('^+z), (1) 

|cunha@mpia.de] 



where T^j'^g is the temperature of the CMB at z = 0, TJ^g = 
2.73 K, sets a fundamental minimum temperature of the 
ISM (assuming local thermal equilibrium, LTE). The in- 
crease of the minimum ISM temperature with redshift af- 
fects the physical conditions of the dust and molecular 
gas in galaxies, boosting both the dust continuum emis- 
sion and_JheJme_Juinm (e.g. [S ilk & Spaans 1997j; 
iBlaiiil [Tggi FCombes et alJ IT999I: iRighi et al.. .2008,) . How- 
ever, at high redshift, the CMB also becomes a stronger 
background against which both the dust continuum and 
line fluxes are measured (e. g. in [C ombes et al. 19991 
Papadopoulos et al. 2000; Obre schkow et al.ll2009t iLidz et all 
|201 1.; Munoz & Furlanetto 20131) . These two competing ef- 
fects have been discussed previously, however only in the 
context of CO line emission (e.g. in ICombes et"al] [19991; 
lObreschkow et al.l 120091 and briefly in other wor ks cited 
above, sometimes incorrectly: Silk & Spaans| [T997b . A full 
exploration of the effects of the CMB on both continuum and 
line emission observations is needed now that ALMA is in 
operation. In this paper, we quantify in detail the implica- 
tions of the CMB on the interpretation of (sub-)millimeter 
observations in terms of galaxy intrinsic dust and gas prop- 
erties. Deriving the correct intrinsic galaxy properties for 
(sub-)mm observations, such as the star formation rate (via 
the total infrared luminosity) and the molecular gas mass 
(via the CO line luminosity or the dust mass) is crucial to 
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our understanding of t he star formation efficiency of galax- 
ies athigh-redshift (e.g. lDaddi et al.l201(]tlGenzel et alJr201(]|: 
iMagdisetalJlMl) . 

In this study, we wish to analyze solely the effect of the 
CMB background at high redshift on (potential) observations 
of 'Milky-Way-like' galaxies at high redshift. Therefore, we 
fix the intrinsic properties of the galaxies under consideration 
and assume no evolution in dust properties and stellar radia- 
tion field, in order to isolate the CMB effects. We note, how- 
ever, that it is likely that high-redshift galaxies have consid- 
erably different properties than low-redshift galaxies of the 
same luminosity. For example, it is likely that galaxies with 
cold ISM such as the Milky Way do not exist due to the harder 
radiation fields in low metallicity environments. Also, the 
physical properties of dust grains may be different at high red- 
shifts. These differences are beyond the scope of this paper. 

In Section m we describe the effects of the CMB on the 
dust continuum of galaxies and the implications on the de- 
rived physical properties of the dust; we also provide a recipe 
to take the effects of the CMB into account when comparing 
models with observations. In Section[3] we analyze the effect 
of the CMB on the observed CO line emission of galaxies, 
focusing on the simplest LTE case, and two more general ex- 
amples of non-LTE cases. Our conclusions are summarized 
in SectionlD 

2. EFFECT OF THE CMB ON (SUB-)MM CONTINUUM EMISSION 

2.1. Dust temperature 

We consider a galaxy at z = 0, where diffuse dust is being 
heated by the radiation field produced by stars in the galaxy, 
and the effects of the CMB radiation are negligible. We as- 
sume that the dust is in thermal equilibrium with the radia- 
tion field, with an equilibrium temperature T ^^^ = 1 8 K, typi' 



calof the diffuse IS M of the Milky Wa y (e.g. fda Cunha et all 
l2008t see also, e.g. iGroves et al.ll2012l for an analysis of the 
cold dust in M31). If we place exactly the same galaxy (i.e. 
same starlight intensity heating the dust, same dust properties 
etc.) at z = 6, where the CMB temperature has increased from 
2.73 K to 19.1 K, i.e. higher than T^^^ (eq.[B, how does this 
higher CMB temperature affect the dust temperature of the 
galaxy? The dust grains will absorb the CMB photons, so the 
temperature of dust in the galaxy at z = 6 will be higher than at 
z = even if the other properties of the galaxy are exactly the 
same. In the following, we derive the temperature of the dust 
in a galaxy at redshift z, 7dust(z), when heating by the CMB is 
taken into account. 

If the dust grains are in thermal equilibrium, then they emit 
energy at the same rate they it absorb it, i.e.: 

dt dt 

For a single dust grain, the energy loss rate through emission 

is: dE^^^, 

=47r / dvBi,[Ti^^i(z)]-Ka^Q^^{v,a), (3) 



(4) 



dt JO 
while the rate of energy absorbed per grain is: 

dt 

where 7dust(z) is the equilibrium temperature of the dust 
grains at redshift z, and 2em('^,fl) and Qabsit^,o) are, respec- 
tively, the emission and absorption coefficient s at frequency 
v, for a gr ain of effective radius a (see, e.g., ISpitzerlll978l: 
|Draine"201 1 for more detail). In general, Qem(i',d)=Qshs(i',a) 
dPraine & Lee.. 1984.) . and we assume that the dust properties 




Figure 1. Variation of the dust temperature (red; eg. 1121 and CMB tempera- 
ture (blue; eq. [T] with redshift. For dust grains in a galaxy with equilibrium 



temperature at z = 0, TF 



: 1 8 K, the effect of additional dust heating by the 



CMB starts being non-negligible around z ~ 4. 



are invariant with redshift, i.e. that 2abs(i^, «) does not change 
with z. 

At any redshift, we can write the intensity of the radiation 
field heating the dust as (e.g. Rowan-Robins on et al..,19 79): 



(5) 



where the first term, J*iz), corresponds to the contribution by 
the radiation produced by stars in the galaxy, and the second 
term, Bi/[7cmb(z)], is the contribution by the CMB radiation, 
which emits as a black body of temperature Tcmb(z)- Thus, 
at a given redshift, using eq. |5] the equation of the energy 
balance of the dust grains in thermal equilibrium (eq. |2) takes 
the form: 

/ diy Q.^hs(i^,a)J*(z)+ diy Q.^hs('^,a)B„[TcMB(z)] 
Jo Jo 



diy Q^bsi'^,a)B„[Tda^tiz)]- (6) 



From this equation, the total energy produced by stars in the 
galaxy (that is absorbed by the dust grains) can be written as: 



diy Q-ibs.(v,a)Jl(z) 



diy Q,bs(y,a){B,m,,,(z)]-BATcMB(z)]}. (7) 

We assume here that the properties of the galaxy do not 
change with redshift, i.e. the stellar radiation field, J*, re- 
mains the same. Therefore, for any redshift z. 



di^Q,Ui^,a)J:(z)= / di^Q,W^,a)j;(z = 0), (8) 
/() Jo 
and combining with eq.|2]we can write: 



dj^ Q,b,W,a){BATdu.i(z)]-BATcMB{z)]} 



diy QaU^y, a) {BAT^:^:]-B,[T^^\,]} 



(9) 



At (sub-)mm wavelengths, the grain emissivity can be 
approximated with a power-law function with frequency. 
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Qabs(i^ifl ) i^^, where /3 is the so-called dust emissivity in- 
dex (e.g. lDraine & Lej|1984l) . Thus, we can replace this term 
in eq.|9]and re-arrange in order to solve for 7dust(z) as follows: 

/•oc poo 

Jo Jo 



z=0 ■ 



The integral of a modified black body v^B^^iJ) is: 



dv vi^B^iT) oc T 



4+0 



(10) 



(11) 



From eqs.[Tl[TO]and[TTl we obtain the following equation for 
the equilibrium dust temperature at redshift z: 



4+/3 



T,udz)={i&'^+(Tizr^[(i+zr^-i\)- . (12) 

In Fig.[Tl we plot the evolution of dust and CMB tempera- 
tures with redshift, assuming a dust temperature of 18 K at 
z = 0. This shows that, for this temperature, the effect of 
dust heating by the CMB becomes non-negligible at z ~ 4; 
the higher the redshift, the dust temperature asymptotically 
approaches the temperature of the CMB. 

2.2. Effect on the intrinsic far-IR/sub-mm dust SEDs 

Fig.[T]shows that the effect of dust heating by the CMB on 
the (cold) dust temperature of galaxies is to first order negligi- 
ble at redshifts z < 4. However, beyond this redshift, the dust 
temperature increases with redshift due to the extra heating by 
the CMB as demonstrated in the previous section. To investi- 
gate how this affects the far-infrared and sub-millimeter SEDs 
of the galaxies, we plot, in Fig.|2] the predicted intrinsic SED 
of dust grains in thermal equilibrium with T^^^l = 1 8 K, at red- 
shifts z= 0.1,2.5,5 and 10. In each panel, the thick grey line is 
identical and shows the un-corrected SED (dust heating only 
by star formation i.e. the radiation field 7*) and the black line 
shows the predicted intrinsic SED when taking into account 
additional dust heating by the CMB. We stress that this is the 
intrinsic SED, not the actually observed one (as discussed in 
Section |23] l. As the dust temperature increases due to increas- 
ing CMB temperature, this affects the dust SEDs in two ways, 
as seen in Fig.|2] (i) the peak of the SED shifts towards lower 
(rest-frame) wavelengths (i.e., higher observed frequencies) 
as the dust temperature increases with redshift; (ii) the total 
luminosity increases by a factor [T^usiiz) /T^^^^]'-'^^^ (eq. fTTTl. 
due to additional energy from CMB photons absorbed by the 
dust grains. 

Finally, to check the impact on the intrinsic (sub-)mm con- 
tinuum flux densities, in the top panel of Fig.|3]we plot the ef- 
fect of the CMB heating on the continuum fluxes in 9 ALMA 
bands from 38 to 870 GHz, as a function of redshift. This 
shows that the CMB has practically no effect on the fluxes 
up to z — 4, but after that redshift it contributes increasingly 
to boost the intrinsic (sub-)mm fluxes in the ALMA bands. 
The effect is strongest for the higher-frequency bands, which 
sample the dust SEDs closer to the peak (see Fig.|2]l. 

2.3. Detectability of dust emission against the CMB 

background 



Extra dust heating as described in Section lZ2] is not the only 
effect of the CMB radiation on the observed (sub-)millimeter 
fluxes. For any (sub-)mm galaxy observation, the dust con- 
tinuum is always measured against the CMB. Here we dis- 
cuss how this affects the detectability of the dust continuum 
of galaxies at high redshifts. 

We assume that the stellar radiation can be neglected at 
(sub-)millimeter wavelengths, i.e. 7* = 0. From solving the 
radiative transfer equation, the (rest-frame) intensity per unit 
frequency from a galaxy at redshift z in the (sub-)millimeter 
is: 

/, = [l-exp(-T,)]B,[rdust(z)] + exp(-T,)B,[rcMB(z)], (13) 

where Ti, is the dust optical depth. The first term of this equa- 
tion corresponds to the emission by dust in thermal equilib- 
rium with temperature rdust(z) (given by eq. [T2] |. and the sec- 
ond term is the contribution from CMB radiation that is trans- 
mitted through the galaxy ISM (i.e. the fraction of the un- 
derlying CMB that is not absorbed by dust). We assume that 
the dust is optically thin in the (sub-)millimetefl i.e. <^\, 
then exp(-rix) ~ I-t^- The optical depth can be expressed 
as: = S^/Ky, where is the surface mass density of dus^^ 
(in units of g cm"^) and is the mass absorption coefficient 
(cross-section per unit mass, in units of g~' cm"^). Therefore, 
eq.[T3]can be re-written as: 



Iv = SrfK^By[rdust(z)] + (1 - ^di^u)BATcMYi(z)] ■ 
The observed flux of the galaxy is: 



70bs 



■ n- 



(14) 



(15) 



(l+z)3 



where il is the solid angle subtended by the galaxy, which 
can be written in terms of the physical area of the galaxy A 
and the luminosity distance c/l, as: $7 = (1 +z)'^A/dl^. From 
eqs.[T4land[T5]we obtain: 

^ u [ 1 

SrfK^B^[rdust(z)]+(l-SrfK^.)B^[rcMB(z)] 

^ (16) 
The flux is always measured relative to the CMB back- 
ground. The flux measured against the CMB is obtained by 
subtracting the observed CMB flux (over the solid angle ^Y) 
from the predicted flux in eq. [16] We obtaiiJ3 

j;,obs against CMB _ j;,obs OD Ft 2=0 1 

^i//(l+z) - ^u/{l+z) - ^^tSuHl+<) L-fCMBj 

_ obs _ o [^CMb(z)] 
-^Wd+z) (1+^)3 



= -^j—^T^dK^ [rdust(z)] -By [rcMB(z)] ) 

= —^MdKyiBy [7dust(z)] -By [7cmb(z)] 



(17) 

Based on this equation, the fraction of the intrinsic dust 
emission from the galaxy (i.e. the first term of eq.fTJt that we 

' This is a reasonable assumption, as the dust column required for the 
galaxy to be optically thick at (sub-)millimeter wavelengths are unrealis- 
tically high. For example, using a dust mas s absorption coefficient of 
Ki/(850 fim) = 0.77 g"' cm^ IDunne et aU2QO0l) . a galaxy of total dust mass 
10** Mq would need to have a radius smaller than 100 pc to be optically thick 
at 850 fJLm. 

^ This can be written in terms of the total dust mass and the physical 
area of the galaxy A: = Mj /A. 

^ The observed CMB intensity is related to the rest-frame CMB intensity 
as: S./(i+,)(r(?^V = 5='[7fcMB(z)]/(l +1?- 
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Figure!. Spectral energy distributions of dust grains in thermal equilibrium with T^^^ = 18 K, with total luminosity powered by starlight in the galaxy of 
10 Lq, at redshifts z = 0.1,2.5,5 and 10 (redshift is indicated in the upper-light comer of each plot), grey line - intiinsic SED when no dust heating by the 
CMB is included; black solid lin e - i ntrinsic SED including dust heating by the CMB; black dashed line - actual measured SED, when contrast with the CMB 
background is included (Section [2.31 . The colored vertical lines indicate the observed frequency (and rest-frame wavelength) sampled by the 9 ALMA bands 
(Uabs = 38, 80, 100, 144, 230, 345, 430, 660 and 870 GHz). 



can actually measure against the CMB at a given frequency 
i^obs = i^/(l +z) is given by: 



rrobs against CMB 
/^intrinsic 

1^/(1+2) 



^ Bu[Tcmb{z)\ 
B^[TaUz)\ ■ 



(18) 



If 7dust(z) = 7cmb(z), i e. the dust is in thermal equilibrium 
with the CMB, then no intrinsic flux can be detected against 
the CMB; on the other hand, if 7dust(z) ^ 7cmb(z) (usually the 
case at low-z), then practically all the intrinsic flux is detected 
against the CMB. In the bottom panel of Fig.[3l we plot the 
fraction of flux observed against the CMB as a function of red- 
shift in the 9 ALMA bands from 38 to 870 GHz, for a galaxy 
with intrinsic dust temperature T^^^l = 18 K. This figure shows 
that, as the redshift increases and the CMB temperature gets 
closer to the dust temperature (as shown in Fig. [1}, it becomes 
increasingly difficult to detect the dust emission against the 
CMB, and this effect is higher for the lower- frequency bands. 
In Fig. 121 we add a black dashed line that represents the actu- 
ally observed SED of a galaxy with dust intrinsic temperature 
-^dust' = 18 K at different redshifts. This is obtained by mul- 



tiplying the intrinsic SED emitted at each redshift (including 
heating by the CMB; black solid lines in Fig.[2]i by the factor 
given by eq.[T8]at each redshift. Figs.|2]and|3](bottom panel) 
show that the fraction of flux from a galaxy at a given redshift 
that can be detected against the CMB varies with frequency. 
This can have severe implications on the interpretation of ob- 
served (sub-)millimeter SEDs at high redshift in terms of dust 
properties. For example, for the SED in Fig. |2] at z = 5, the 
Rayleigh-Jeans tail of the dust emission observed against the 
CMB (black dashed line) is much steeper than in the intrin- 
sically emitted SED (black solid line), which would be inter- 
preted as a higher emissivity index /3 of the dust grains, and 
would change the inferred dust luminosity, mass and temper- 
ature, as we discuss in more detail in the next section. 

2.4. Discussion: effect on the physical interpretation of the 
observed (sub- )mm dust emission 

In the previous sections we show that the CMB has an im- 
pact on the peak of the observed dust SED and its slope in the 
Rayleigh-Jeans regime. Here we investigate how this may 
affect the inferred dust temperature, luminosity, emissivity 
index and mass from typical (sub-)millimeter observations. 
It is usual practice to interpret observed far-infrared/sub-mm 
galaxy SEDs by comparing them with modified black body 
(MBB) functions (of the form i''^ B^[Tdust]) that describe the 
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Effect of extra dust heating by the CMB 
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Figure 3. Top panel: Effect of additional dust heating by the increasing CMB 
temperature with redshift. We plot the ratio between the predicted flux den- 
sity in each ALMA band when including dust heating by the CMB and not 
including dust heating by the CMB, for dust grains with an equilibrium tem- 
perature of 1 8 K at z = 0. This figure does not include contrast with the CMB 
background. This shows, as in Fig. |2] that the highest frequency bands ai'e 
the most affected by dust heating by the CMB. Bottom panel: Effect of the 
CMB as an observing background. We plot the ratio between the flux in 
each ALMA band that can be measured against the CMB at a given redshift, 
and the intrinsic flux emitted by the galaxy at that frequency. As for the top 
panel, we assume T^^^^ = 1 8 K and, at each redshift, include the extra heating 
contributed by the tMB in the intrinsic emitted flux. 

emission by dust in thermal equihbrium using only three pa- 
rameters: the dust e missivity index /3, temperature rdnsf, and 
total luminosity (e.g. lDunne et alj |200Q';'Kla as et alJlioOlh . 

We consider the fiducial model discussed in the previous 
sections, with an intrinsic dust temperature of T^^^^ = 18 K, 
dust emissivity index (3 = 2, total luminosity 10" L©, and total 
dust mass 1 .6 x 10^ Mq. We assume that our model galaxy is 
at z = 5, and compute the predicted fluxes in five typical (sub- 
)millimeter bands at 100 GHz (3 mm), 230 GHz (1.3 mm), 
300 GHz (1 mm), 345 GHz (870 Aim) and 670 GHz (450 ^m). 
These fluxes are supposed to be the 'real' observed fluxes, 
therefore we take into account the effect of extra heating and 
extra background provided by the CMB at z = 5, as discussed 
in Sections [2. ll to l2.3l We assume a typical flux uncertainty of 
5% in each band. In order to test how biased our temperature, 
emissivity, luminosity and dust mass estimates for this model 
galaxy would be if we ignored the effect of the CMB, we fit 
the model fluxes in the five (sub-)millimeter bands using a set 
of models that does not include the effect of the CMB; we also 
compare with a set of models that correctly includes the effect 
of the CMB. In practice, we compare our synthetic observed 



fluxes in five bands to a grid of modified black bodies where 
we vary the intrinsic dust temperature T^^^^ between 10 and 
30 K, and the emissivity index (3 between 1 and 3. 

In Fig. m we plot the results of this fitting when includ- 
ing the CMB effects (in black) and the results when ignor- 
ing the additional effects of the CMB (in red); for each case, 
we plot the best-fit SED and the marginalized probability 
distribution functions (PDFs; com puted using the approach 
described in Ida Cunha et al.l 120081) for the dust temperature, 
emissivity, total luminosity and mass, with the best-fit value 
for each parameter (i.e. the value that minimizes y^) indi- 
cated by a vertical dashed line. We first check that the best- 
fit model parameters obtained when including the CMB ef- 
fect (indicated by the black dotted lines) are equal to the in- 
put parameters - this confirms the robustness of the method. 
When the SED is fitted with models that do not include the 
effects of the CMB, there is a significant difference between 
the best-fit model parameters and the input (i.e. 'real') pa- 
rameters. The intrinsic dust temperature and emissivity index 
are overestimated {T^^^^ = 19 K instead of 18 K, and (3 = 2.1 
instead of 2.0), while the dust luminosity and the dust mass 
are underestimated (log(L(just/i'0) = 10.87 instead of 11.00, 
and log(Mdust/A^0) = 7.35 instead of 8.20). The effect on the 
estimated dust mass is the strongest because the dust mass de- 
pends on the luminosity, temperature and emissivity index as 
Mdust oc Ld ust T7 \'^^^\ This confirms what was already hinted 
in Section 12.31 the CMB makes the SEDs look hotter and 
steeper, and Fig.|4]shows that this has a significant impact in 
the deduced properties of dust in high-redshift galaxies when 
fitting standard MBB models to the observed dust emission. 
A steeper emissivity index would imply different dust proper- 
ties at high redshift, and the difference in inferred dust masses 
would also change our understanding of ISM enrichment by 
dust in high-redshift galaxies. 

In Fig.|5] we investigate these effects in a more systematic 
way. We use the same model and perform the similar modified 
black body fits as in the test described in Fig.Hl but at various 
redshifts between z = and z = 10, and for two intrinsic input 
dust temperatures, T^^f = 18 K (in blue) and 7^;=° = 40 K (in 
red). The four panels of Fig. |5]show the difference between 
the best-fit parameter value when ignoring the CMB effects 
on the SEDs and when including the CMB effects, for the 
dust temperature (a), the dust emissivity index (b), the total 
luminosity (c), and the dust mass (d). We show, for each tem- 
perature, the result of fitting the same five (sub-)mm bands 
considered in Fig. |4l observed frequencies i^'obs = 100, 230, 
300, 345 and 670 GHz. Not surprisingly, the overestimation 
of temperature and emissivity index and the underestimation 
of luminosity and mass are worse at high redshifts, and for the 
lowest intrinsic dust temperature, since in this case the tem- 
perature contrast between the CMB and the dust emission is 
small. The cold dust mass can be severely underestimated, 
up to about two orders of magnitude at z — 10. We note that, 
at very low redshifts, we underestimate the dust temperature 
and luminosity, and overestimate the emissivity index, even 
though we expect the CMB to have a minimal effect at low 
redshifts: the differences should be close to zero in all panels 
at z = 0. This offset is explained by the fact that, when fitting 
fluxes at observed-frame frequencies 670 GHz and lower, as is 
the case for the solid lines in these plots, we sample the SED 
significantly lower frequencies from its peak, and relatively 
close to the peak of the CMB emission at low redshifts. 
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Figure 4. Comparison between fits to our fiducial model SED (intrinsic dust temperature of ^ = 18 K, dust emissivity index /3 = 2, total luminosity lO" Lq, 
and total dust mass 1.6 X 10^ Mq) at z = 5. when including the effect of the CMB in the fit (in Black) and when ignoring the effect of the CMB in the fit (in red). 
The 'observed' fluxes are indicated as blue diamonds in the top-left panel; we use five common (sub-)millimeter bands at 100, 230, 300, 345 and 670 GHz. The 
grey SED is the intrinsic SED, i.e. before CMB corrections - the actual observed SED i.e. with the effects of the CMB included as for Fig. |2] is plotted as a 
dashed black line. The best-fit SED when fitting the blue fluxes but ignoring the effects of the CMB is plotted in red. The top-right panels and the two bottom 
panels show the probabil ity density functions (PDFs) of the constrained parameters (dust temperature, emissivity index, luminosity and mass) computed using a 
Bayesian fitting method 4da Cunha et aU2008ft . with the best-fit value for each parameter indicated by a vertical dashed line. 



2.5. How to account for the CMB when interpreting real 
continuum measurements 

We have shown in the previous section that, when inter- 
preting real (sub-)nim continuum measurements in terms of 
intrinsic dust properties, we must take the effect of the CMB 
into account. This is the case for comparing observations with 
dust emission models, such as modified black bodies (as de- 
scri bed above), or even more complex dust emission models 
(e.g. Ida Cunha e t al. 2008; Draine & Li 2007). To summarize, 
we can estimate the intrinsic properties of the dust in an ob- 
served galaxy at redshift z by following these steps: 

1. for a given model with intrinsic temperature Tl^^^ and 
emissivity index /3, compute the extra heating provided 
by the CMB using eq. [121 the actual dust temperature 
is then rdust(z); 

2. include the extra luminosity provided by the CMB 
heating, i.e. multiply the dust emission model by 
[rdus.(z)/7d;=°](4^«; 

3. compute dust emission model in the observed frame, 

pintiinsic . 

4. to account for the effect of the CMB as an observ- 
mg background, obtam ^^.y^j^^j by multiplying 
the dust emission model by the factor given in eq. [Tsl 

i-{B.[rcMB(z)]/B.[rdust(z)]}; 

5. compare the modified dust emission model (parameter- 
ized in terms of the intrinsic dust properties from step 
1) directly with the observations. 



3. EFFECT OF THE CMB ON THE CO LINE EMISSION 

We now discuss the effect of the CMB on the CO emis- 
sion at high redshifts. The mechanisms by which the CMB 
affects the CO excitation are not as straightforward as for the 
dust case discussed in the previous section, since the local 
thermal equilibrium (LTE) conditions can only be assumed in 
particular cases, and in reality molecular clouds are often in 
non-LTE conditions. In general, the effect of the CMB on the 
detectability of high-redshift CO lines is two-fold, as in the 
case of dust emission. On one hand, the higher temperature of 
the CMB at higher redshifts can affect the excitation of CO by 
helping populate high rotational levels (through mechanisms 
that we discuss below), thus increasing the line luminosi- 
ties fo r high-rotational number transitions (e.g. ISilk & SpaansI 
Il997h . On the other hand, as the CMB temperature increases, 
it becomes a more important observing background a gainst 
which the CO lines must be detected (see also Combe s et al.l 
[T999tlPapadopoulos et al.ll2000l:IObreschkow et al...2009lK 

For simplicity we assume throughout this section that the 
gas and dust are efficiently coupled, such that the dust temper- 
ature, Tdust, sets the minimum kinetic temperature of the gas, 
Tkin- This implies that the gas heating via collisions with dust 
grains is 100% effective (Tielens 2005). This can be the case 
for high optical depths (Tielens & Holle nbachll 19851) . but, in 
most realistic cases, the heating of gas by dust may not be 
100% effective. However, even in less optimal conditions the 
gas and dust temperatures are coupled so that if the dust tem- 
perature increase s, the gas temperature increases by a similar 
amount (Tielensi 120051) . We discuss deviations from the as- 
sumption that Tkin = Tdust where relevant. 

3.1. Some general aspects of CO excitation 
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Figure 5. Difference between different dust parameter estimates when not including the CMB effects minus when including the CMB in the fits (fitting the fluxes 
in five bands with fobs = 100,230,300,345,670 GHz), for the dust emission model with /3 = 2 and intrinsic temperature = 1 8 K (in blue) and rj~~i =40K(in 
red): (a) temperature; (b) emissivity index; (c) logarithm of the total luminosity; (d) logarithm of the total mass. This shows that, as in Fig.|4] the dust temperature 
and emissivity index tend to be overestimated when ignoring the effect of the CMB on the SEDs at high redshift, and the total dust luminosity and mass are 
underestimated. As expected, these systematic effects are less severe for higher intrinsic dust temperatures. 



Here we summarize some of the main asp ects of CO ex- 
cit ation; more detai l s can be found e.g. in ISpitzerl (11978)) 
or Ivan der Tak et al.l (120071) . The frequency of the photon 
emitted during the rotational transition of the CO molecule 
from a level y„ to a lower-excitation level Ji = Ju- I is: 
i^„i = [E{Ju)-E(Ji)]/h = hJii/iAn^mr^) ~ i^coJu, where/; is the 
Planck constant, m is the reduced mass of the CO molecule, 

is the effective distance between the two atoms, and i^co = 
115.2712 GH2Q. 

The transition between rotational energy levels of CO 
in molecular clouds occurs mainly via collisions with H2 
molecules, with rates given by the collision coefficients C,,/ 
(for collisional de-excitation) and C/„ for (collisional ex- 
citation). These collision coefficients depend mainly on 
the velocity-integrated collision cross-section (which depends 
among other things on the kinetic temperature of the gas, Tkin), 
and the number density of collision partners, in this case hh, ■ 
Photons of frequency i^ui are emitted when a molecule goes 
from an upper energy level 7„ level to a lower energy level 
Ji via spontaneous emission, which occurs at a rate given by 
the Einstein coefficient A,,/, or by stimulated emission, which 
occurs at a rate B„/t7, where B„/ is the Einstein stimulated 
emission coefficient and U is the line profile-weighted mean 
energy density of the radiation field. Photons of the same fre- 

■* In reality, the frequency of /„ > 1 transitions are close but not exactly 
multiples of Uco because rotation of the molecule makes re vary slightly (due 
to centrifugal forces); the exact frequencies of the transitions can be found in 
[http ;//splatalogue.net 



quency from the radiation field are absorbed at a rate B/„t7, 
where B/„ is the Einstein absorption coefficient. The Einstein 
probability coefficients are related by: 



giBiu = guBui = 



giAui 



(19) 



where gu and gi are the statistical weights of the levels (given 
by 2Jii + 1 and 27/ + 1 = 27,, - 1, respectively). The excitation 
state and line absorption/emission by CO molecules is com- 
puted by taking into account collisions and radiative transi- 
tions. In statistical equilibrium, the populations of levels u 
and I are given by: 



n„(A„, + B„,t/ + Cui) = ni(B,M + Ciu) , 



(20) 



where U = 4ttI/ c is the average energy density of the radiation 
field, and / is the intensity of the radiation field 7^ averaged 
over the line profile. We can assume I^, = 8,^(7^^^), where TJad 
is the radiation temperature. 

The 'excitation temperature' of a transition from 7„ to 7/, 
^ex'c' by definition, characterizes the level populations result- 
ing from eq.|20]using the Boltzmann equation: 

riu gu [ \ .Tlx 

= —ex^y-j^] . (21) 



ni gi 



kT-'" 



If collisions dominate the excitation of CO, then the level 
populations are set by the collision coefficients and the kinetic 
temperature, i.e. 

Qu "h gu ( hv \ 

7^ = — = — expl-T— 1. (22) 

Cui ni gi V kTi^n'' 
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The radiative transfer equation is: 

— =-Kj,y + ey, (23) 
as 

where = (hv / c){niBiu-nuBui) is the absorption coefficient, 
and ei, = {hv / ATi)nuAui is the emissivity. The optical depth of 
a transition is defined as: 



and we can define the source function Su as: 

€^ c yiuAiii 



4tt riiBiu-riuBiii 
and re- write the radiative transfer equation (eq. E 

which has the solution: 

= [1 - exp ( - T^')] 5, + /° exp ( - tI;'\ 



(24) 
(25) 
(26) 
(27) 



The first term is the 'intrinsic' line flux, and the second term 
is the transmitted background intensity. From eqs. |20]|2n and 
|25]follows that: 

S.= ^ = B, (7;4) . (28) 



exp 



hu 



■1 



If we assume the only background is the CMB, then the (rest- 
frame) intensity is: 

U=[l-exp{-Ti-')]B,{T^^;)+exp{-Tl-)B,[TcMB(z)] , 

(29) 

and, similarly to eq. [15] the observed velocity-integrated flux 
of the line is: q 

S»^^l=^^I^. (30) 



The velocity-integrated flux of the line measured against the 
CMB background is then (see also eq.|17|i: 



[obs against CMBJ 
■/(1+z) ' 



[obs] 



3 [l-exp( 



(31) 
that is 



The fraction of intrinsic line flux (first term in eq. 
observed against the CMB is: 



cJi, [obs against CMB] 
r-/,, [intrinsic] 



B,[Tcmb(z)] 



(32) 



This is similar to eq. [18] and reflects the fact that we cannot 
detect the emission from lines where the excitation temper- 
ature is the sa me as the background rad i ation (as also in- 
cluded in. e.g. [Scoville & Solomon|[i974t Ivan der Tak et al.l 
120071 lObreschkow et al.ll2009l]^ 

The previous equations show that the observed intensity of 
a CO line depends on the optical depth t^" and the excitation 
temperature Tj^^ of each transition, which are set by the sta- 
tistical equilibrium equation (eq. |20] | and the radiative trans- 
fer equation (eq. |26] l. Computing t^" and T^^^ is a non-trivial 
problem since the level populations and the radiation field are 

^ We note that CO is unlilcely to be obsei'ved in absorption against the 
CMB at any redshift, as the ubiquitous CMB radiation will effectively heat 
any gas (either via dust or line heating) to the temperature of the CMB. 



coupled. A commonly used method to deal with this is the 'es- 
cape probability method' which gives the probability (3 that a 
photon emitted in the transition escapes from the cloud. In 
this case, the radiation field becomes 0(1- (3), and the way 
f3 depends on the optical depth tI" is specified by the geome- 
try. The most widely used approximation is the large velocity 
gradient (LVG) approach, which assumes radially expanding 
spheres characterized by a velocity gradient dv/dr along the 
line of sight ( e.g. Sco ville & Solomonll974tlWeiB et al.ll200l 
Ivan der Tak "et al.ii2007lK 

3.2. The LTE case 

If collisions dominate, i.e. C,,/ ^ A^i (this is the case for 
high densities), then Tj^^. = Tj^in for all transitions. This is the 
local thermal equilibrium (LTE) case, and it implies that all 
the lines are thermalized (e.g. Obreschkow et al. 2009). As 
mentioned above, we assume thermodynamic equilibrium be- 
tween the gas and the dust, i.e. Tkin = Tdust- Since in LTE, 
^ex'c = ^kin, then TJ^^ = Tdust, and therefore the excitation tem- 
perature increases with redshift in the same way as the dust 
temperature (eq.[T2]i. 

In Fig. [6] we plot the ratio between the line luminosity ob- 
served against the CMB and the intrinsic luminosity for the 
CO transitions from 7„ = 1 to y„ = 10 using eq.[32] We adopt 
two different intrinsic gas kinetic temperatures, T^Q^^ = 18 K 
(top panel) and Tj^J^" = 40 K (bottom panel). This figure 
shows that the decrease of CO line flux due to the CMB back- 
ground can be a very significant effect, substantially decreas- 
ing the measured line flux for sources at high redshifts (see 
also lCombes"etal. 1999: Obreschkow et al. 2009 ). The effect 
of the CMB background is more pronounced for colder in- 
trinsic excitation temperatures (which are closer to the CMB 
temperature at any redshift). This has strong implications for 
the detectability of CO lines from 'quiescent' galaxies at high 
redshift, where the radiation field from star formation in the 
galaxy may lead to relatively low Tidn (and therefore Tj^^.). For 
TjJ^" = 18 K, typically less than half of the total line luminos- 
ity can be measured at z > 4, with even lower fractions for the 
lowest-/,, transitions. For higher gas excitation temperatures, 
the effect is not as strong. For 7j?J^" = 40 K, over 70% of the 
intrinsic flux is recovered. The decrease of measured CO line 
flux at high redshift shown in these plots has strong impli- 
cations for measuring the molecular gas (Mh,) mass of high- 
redshift galaxies usi ng the CO line luminosity , Leo- Since 
Mh, oc Leo (e.g. Solo mon & Vanden Boutll2005l) . if we do not 
correct for the CMB in order to obtain the intrinsic CO lu- 
minosity (using eq. [32] ), then we may severely underestimate 
Leo and consequently underestimate the molecular gas mass 
of the galaxies (by more than a factor of 2 at z > 4 for 18 K 
gas). 

Fig. [6] also shows that the fraction of CO line luminosity 
that is measured against the CMB for a given gas excitation 
temperature is not the same for all CO lines as a function of 
redshift. This implies that the shape of the observed CO spec- 
tral line energy distributions (SLEDs) will change with red- 
shift, even if the intrinsic shape, i.e. emitted SLED (which 
depends only on the gas properties and excitation) does not 
change. Like in the case of the dust SEDs, this can affect how 
one interprets the properties of high-redshift galaxies based 
on their observed CO SLEDs, namely the excitation condi- 
tions of the gas. We investigate how large this effect is by an- 
alyzing the CO SLEDs is as a function of intrinsic excitation 
temperature and redshift. We consider three intrinsic kinetic 
temperatures, T^f = T^^^ = 18, 40, and 100 K. For each tem- 
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Figure 6. Ratio between the line (velocity-integrated) flux observed against 
the CMB background and the intrinsic line flux for different CO transitions 
from J„ = 1 to /„ = 10 (shown with different colors). In the top panel, the 



intrinsic gas kinematic temperature is TJ 



18 K; in the bottom panel, we 



show a case with higher intrinsic gas temperature, T^^ = 40 K. Since this 
is the LTE case, the excitation temperature Texc is the same for all levels 
and equal to Ti^i,,. These plots include both effects of increasing CMB tem- 
perature with redshift on observed CO lines: (i) the increase of gas kine- 
matic temperature (and hence the excitation temperature) as described by 
eg. 1121 (ii) the growing importance of the CMB as an observing background 
with redshift. When SJ' [observed against the CMBJ/SJ' [intrinsic] = 1, the 
measured line flux is the total intrinsic flux emitted by the source; when 
5;i' [observed against the CMB]/5^' [intrinsic] = 0, the hne cannot be distin- 
guished from the CMB background. 



perature, we build the expected CO SLEDs by computing the 
velocity-integrated flux of each 7„ > 1 transition normalized 
to the velocity-integrated flux of the 7„ = 1 transition at z = 
(this normalization highlights the fact that the lines become 
less bright due to lack of contrast with the CMB as shown in 
Fig.|6]). We use the first term of eq.|29]to compute the intrin- 
sic CO line luminosity of each transition, with the excitation 
temperature 7 l{v = Tkin, and the l i ne opt ical depth as given by 
equation 5 in lObreschkow et alj (l2009h (derived for the LTE 
approximation): 



T;;"[LTE] = 7.2Tcexp 



2 k Tg/c 



sinh 



2 k Tgx'c 



(33) 



where is a constant. Following Obreschkow et al. ( 200^, 
we fix Tc = 2 (which they show provides a good fit to observed 
CO SLEDs). We compute the expected CO line SLEDs for 
three intrinsic kinetic temperatures at different redshif ts by in- 
cluding both the extra heating by the CMB that increases Tdust 




Figure 7. Variation of observed CO SLEDs with redshift for three different 
intrinsic kinetic temperatures: T^^ = 18 K (blue), 40 K (red), and 100 K 
(yellow). For each intrinsic temperature, we consider four different redshifts 
of the emitting galaxy: 7, = (solid lines); z = 2 (dashed lines); z = 5 (dot- 
dashed lines); and z = 10 (dotted Unes). For each temperature, all the CO 
SLEDs are normalized to the flux of the CO(l-O) line at z = 0, S^iz = 0). 



and hence 7i;i„ and T^J;,, and the brighter CMB background 
as discussed above. The results are plotted in Fig. |7] which 
shows two main effects: (i) the overall flux of the lines de- 
creases at higher redshifts, as discussed above; (ii) the de- 
crease in line intensity is not the same for all transitions, 
which slightly changes the shape of the SLEDs. At fixed 
^kln*'' there is a slight shift of the CO SLED peak towards 
higher J„ with increased redshift. This is a consequence both 
the slightly higher excitation temperature (due to additional 
CMB heating) and the fact that the CMB background affects 
the low-7„ line measurements more than the high-y„ lines (see 
Fig.|6l). As expected, these effects are stronger for the lowest 
Tjli^"; gas with = 18 K is practically undetectable atz > 5. 



3.3. Non-LTE examples 

In real galaxies, the molecular gas is not likely to be in LTE, 
and so the radiative transfer of the lines and the population 
leve ls ha ve to be solved simultaneously as described in Sec- 
tion [TT] In this section, we show how the CMB affects the 
detectability of the CO lines, where the line fluxes and opti- 
cal depths are computed using the LVG model of iWeiB et al.l 
(2005), which assumes an expanding sphere geometry with 
velocity gradient dv/dr. The main free parameters of this 
model are: 

• The kinetic temperature of the gas, Tkin. As in the previ- 
ous section, we assume 7kin(z) = 7dust(z), i-c. that there 
is a 100% effective coupling of the dust and gas tem- 
peratures. We analyze three values of T^^^: 18, 40 and 
100 K. 

• The number density of H2 molecules, n-n^ . We analyze 
two examples: 

- 'low-density' case: n-a^ = 10^^ cm"^ 

- 'high-density' case: = lO"*'^ cm~^ 

• The number density of CO molecules. We assume 
nco = 8 X IQ~^ n^^. 
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Figure 8. Excitation temperatures of the various CO transitions computed using the LVG approximation (using the models of lWeiB eTal . 2005), for a low-density 



model (left-hand panel) and a high-density model (right-hand panel). We consider three intrinsic kinetic temperatures at z = 0: T. 



100 K (yellow). For each model, the kinetic temperature increases with redshift following the dust temperature (eg. 1121 which mcludes extra dust heating by 
the CMB at high redshift). The resulting excitation temperatures at different redshifts are plotted with different line styles: z = (solid), z = 2.5 (dashed), z = 5 
(dot-dashed), and z = W (dotted). We plot only transitions where the optical depth tJ' is higher than 0.1 (Fig.|9). 



18 (blue), 40 (red), and 




• The velocity gradient dv/dr. Assuming that the clouds 
are virialized, we obtain dv/dr=% km s~' pc~' for the 
low density case and dv/dr = 3.9 km s~' pc~' for the 
high density case. 

Using these parameters, the excitation temperature and opti- 
cal depth of each line are computed by coupling the sta tisti- 
cal equihbrium a nd radiative transfer equations (Section 13.11 
IWeiB et al.l l2005). In Figs.[8]and|9l we plot the variation of 
excitation temperatures T^^^ and line optical depths tI", re- 
spectively, with redshift and input kinetic temperature 
resulting from these calculations. These plots show that the 
high-density case approaches LTE for low J„: Tj^^ ~ Tkin up 
to y„ ~ 5 i.e. the low-y„ lines are thermalized (i.e. the H2 
density is higher than the critical density for these levels). For 
higher 7„, Tj^^. < T^in because at these densities and kinetic 
temperatures these levels are not as populated (as they would 
be in the LTE case). We also note that, for the lowest T^^'^, 
at high redshift, the kinetic temperature becomes very close 
to the CMB temperature. Due to the strength of the CMB at 



high redshifts, the CMB defines the minimum excitation tem- 
perature as the CO levels become radiatively dominated. In 
this case, the radiative processes dominate and so the excita- 
tion temperature gets close to Trad i.e. Tcmb- Finally, these 
plots also show that Tj^^ < Tkin at low density and high 7,,, i.e. 
less molecules are in excited states because there is less col- 
lisional excitation. At higher redshifts, 7^^;, is higher because 
Tkin is larger due to the dust being hotter (thanks to extra CMB 
heating) and, at the same time, CMB radiative excitation be- 
coming more important. 

Once the excitation temperature and optical depth of each 
transition are determined using the LVG models, we can com- 
pute the velocity-averaged flux of each line as described in 
the previous sections. In Fig. [TO] we show the ratio between 
the line luminosity observed against the CMB background 
and the intrinsic line luminosity for different CO transitions 
for the low- and high-density case, for intrinsic gas kinetic 
temperatures Tj^Jj"" of 18 and 40 K. These plots can be di- 
rectly compared with the LTE case shown in Fig. |6] As for 
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Figure 10. Ratio between the line lumi nosit y observed against the CMB background and the intrinsic line luminosity for different CO transitions from 7„ in the 
non-LTE examples described in Section [33] (shown with different colors as in Fig|6). In the top panels we show models with intrinsic gas kinematic temperature 



18 K; in the bottom panel, we show the cases with higher intrinsic gas temperature, = 40 K. The excitation temperature for each transition is computed 
using the .Weifi et al.. 1.2005.) LVG model, with the parameters described in Section l33] the left-hand panels correspond to the low-density case, and the right-hand 
panels coiTespond to the high-density case. 




Figure 11. Variation of observed CO SLEDs with redshift for three different intrinsic kinetic temperatures: TJi^ = 18 K (blue), 40 K (red), and 100 K (yellow), 
computed using the WeiB et al. (2005) LVG model for the low-density case (left-hand panel) and the high-density case (right-hand case). For each intrinsic 
temperature, we consider four different redshifts of the emitting galaxy: z = (solid lines); z = 2 (dashed lines); z = 5 (dot-dashed lines); and z = 10 (dotted lines). 
For each temperature, all the CO SLEDs are normalized to the flux of the CO(l-O) line at z = 0, S^iz = 0). 



12 



E. da Cunha et al. 



the LTE case, the contrast between the lines and the CMB 
decreases with redshift, and it decreases more rapidly for 
the low-temperature case, since the temperature is closer to 
the CMB temperature. The main difference between these 
plots and the LTE case are the high-/,, transitions. Since the 
higher-/,, levels are underpopulated compared to the LTE case 
(the lines are not thermalized), their excitation temperature is 
lower than Ikin(z) (see Fig. |8]l, and so the contrast between 
these lines and the CMB is even lower than for the low-7„ lines 
(contrary to what happens in the LTE case). The predicted de- 
crease in contrast is stronger for T^^'^ = 1 8 K, specially for low 
H2 density: in this case, lines with /„ > 5 become practically 
undetectable against the CMB from z = 4. For the lowest-/,, 
transitions, the contrast behaves similar to the LTE case, be- 
cause the low-/,, are thermalized. 

In Fig. [m we plot the CO SLEDs computed using the 
LVG code, i.e. with the excitation temperatures and line op- 
tical depths shown in Figs. [8] and |9] As for the LTE case 
(Fig. IT), these plots show how the flux of the lines decreases 
with redshift due to decreasing contrast with the CMB (and 
also to some extent different line optical depths; Fig. |9]l. In 
this case it is harder to distinguish differences between the 
CO SLED shapes at different redshifts caused by the dif- 
ference in contrast against the CMB or by different intrin- 
sic values of T^^^ and t^" at different redshifts (set by in- 
creasing increasing excitation by the CMB both via extra 
dust heating and radiative processes; Figs. |8] and |9]l, which 
would produce different intrinsic CO SLEDs at different red- 
shifts. However, we have shown that the effect of the CMB 
background is often non-negligible (Figs. [Tol i and must be 
taken into account w hen interpreting observed SLEDs us- 
ing L VG models (e.g. IScoville & Solo monlll974t I WeiB et all 
l200llvan der Tak et al.ll200'7|[ 

We note that if we do not assume a perfect coupling be- 
tween the dust and gas temperatures, i.e. that 7kin(z)=7dust(z), 
the effects presented in this section will be even more pro- 
nounced, i.e. in general the lines look weaker Indeed, if we 
assume that Tjdn to be constant with redshift (i.e. 7kin(z) = 
^kin")' ^'^ h^S^ redshifts we obtain 7kin(z) < 7dust(z), since 
Tdust increases with z due to the extra CMB heating (eq. fTST i. 
In this case, using our LVG modelling, we obtain the result 
that TJ^^ tends to be lower in general (because since r^in is 
smaller the collisions are not as effective bringing molecules 
to excited levels), and so the contrast between the lines and 
the CMB is weaker 

4. SUMMARY & CONCLUSION 

In this paper, we have analyzed in detail how the CMB af- 
fects (sub-)millimeter observations of both the dust contin- 
uum and CO lines in high-redshift galaxies. We have shown 
that, at high redshifts, the CMB becomes an increasingly sig- 
nificant additional heating source for the dust and gas, boost- 
ing their temperatures and enhancing the dust continuum and 
CO line emission of the galaxies. However, we show that 
at higher redshifts, since the CMB is hotter and therefore 
brighter, the contrast of the intrinsic dust (and line) emission 
against the CMB decreases. We have quantified how this af- 
fects the detectability of dust and lines in specific cases with 
fixed temperatures. Additionally, we provide general correc- 
tion factors to compute what fraction of the intrinsic dust (and 
line) emission can be detected against the CMB as a function 
of frequency, redshift and temperature, and a recipe to include 
these effects when comparing dust emission models with ob- 
servations. 



We show that neglecting the CMB effect on (sub-)mm ob- 
servations of high-redshift galaxies can lead to significant dif- 
ferences in the derived intrinsic properties of dust, in partic- 
ular a severe underestimation of the dust mass. Specifically, 
the intrinsic dust mass of cool dust (Tdust < 20 K) at z > 5 can 
be underestimated by at least one order of magnitude. The in- 
ferred total dust luminosity (and hence star formation rate) 
is less drastically affected: between 30% and 50%. Simi- 
larly, if the effect of the CMB as an observing background 
is neglected when measuring CO lines, this can lead to wrong 
interpretations of the molecular gas properties (such as total 
mass, density and kinetic temperature) via the observed CO 
luminosity and the SLEDs. For low gas kinetic temperatures 
(18 K) and different density scenarios, we find that, at z > 5, 
less than 20% of the intrinsic CO(l-O) or CO(2-l) fine fluxes 
can be measured against the CMB. This implies that, without 
proper corrections, the inferred molecular gas mass would be 
20% of the intrinsic value. At higher gas kinetic temperature 
(40 K), the underestimation of the gas mass would be less dra- 
matic since between 20% and 60% of the intrinsic line fluxes 
are measured (depending on the density), however this is still 
significant when studying the gas reservoir/star formation ef- 
ficiency of high-redshift galaxies. 

Finally, our results imply that the cold ISM of galaxies at 
high redshift, with intrinsic temperatures of about 20 K, if 
existent, will be difficult to measure even with the unprece- 
dented sensitivity of modern (sub-)millimeter observatories 
such as ALMA, simply because the contrast between the con- 
tinuum (and line) emission against the CMB background de- 
creases dramaticaUy at z > 4. 
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